/  v 


NPS55-87-012 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


ESTIMATION  OF  THE  PROBABILITY  OF  A  LONG  TIME 
TO  THE  FIRST  ENTRANCE  TO  A  STATE  IN  A 
SEMI-MARKOV  MODEL 


PATRICIA  A.  JACOBS 


L 


SEPTEMBER  1987 


FedDocs 

D    208.14/2 

NPS-55-87-012 


Approved  for  public  release;  distribution  unlimited. 

Prepared  for: 

Office  of  Naval  Research 

Arlington,  VA  22217 


NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFORNIA 


Rear  Admiral  R.  C.  Austin  K.  T.  Marshall 

Superintendent  Acting  Provost 

Reproduction  of  all  or  part  if  this  report  is  authorized. 
This  report  was  prepared  by: 


I  JRITY  CLASSIFICATION  OF  THIS   PAGE 


301 


REPORT  DOCUMENTATION  PAGE 


REPORT  SECURITY  CLASSIFICATION 


NCLASSIFIED 


lb    RESTRICTIVE   MARKINGS 


J  "SECURITY  CLASSIFICATION  AUTHORITY 


DECLASSIFICATION /DOWNGRADING  SCHEDULE 


3     DISTRIBUTION  'AVAILABILITY  OF   REPORT 

Approved   for  public   release;   distribution 
unl imited. 


Performing  organization  report  number(S) 
>S55-87-012 


5    MONITORING  ORGANIZATION   REPORT   NUMBER(S) 


NAME  OF  PERFORMING  ORGANIZATION 

ival    Postgraduate  School 


6b    OFFICE  SYMBOL 
(If  applicable) 

Code  55 


7a    NAME  OF   MONITORING  ORGANIZATION 


ADDRESS  (City.  State,  and  ZIP  Code) 

nterey,   CA  93943-5000 


7b    ADDRESS  (City.  State,  and  ZIP  Code) 


NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

ffice  of  Naval   Research 


ib    OFFICE  SYMBOL 
(If  applicable) 

Code  430 


9    PROCUREMENT  INSTRUMENT  IDENTIFICATION   NUMBER 

N0001487WR24005 


ADDRESS  (City,  State,  and  ZIP  Code) 


10    SOURCE  OF   FUNDING   NUMBERS 


rlington,   VA     22217 


PROGRAM 
ELEMENT  NO 

61153N 


PROJECT 
NO 


RR014-05-01 


TASK 
NO 


4114531-02 


WORK    UNIT 
ACCESSION    NO 


TITLE  (Include  Security  Classification) 

[STIMATI0N   OF  THE   PROBABILITY  OF  A  LONG  TIME   TO   THE   FIRST   ENTRANCE  TO  A  STATE    IN  A  SEMI 
4ARK0V  MODEL 


!.  PERSONAL  AUTHOR(S) 

JACOBS,    PATRICIA  A. 


}a.  TYPE  OF  REPORT 

Fechnical 


13b    TiME  COVERED 
FROM  TO 


14    DATE  OF  REPORT   ( Year.  Month,  Day) 

Sep  1987 


15    PAGE  COUNT 

27 


5.  SUPPLEMENTARY  NOTATION 


COSATI  CODES 


FIELD 


GROUP 


SUB-GROUP 


18    SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Estimation  of  the  survivor  function  of  a   first   passage   time; 
semi-Markov   process;   empirical   moment  generating   function; 
asymptotic  exponential ity;    strong  law;   central    limit  theorem; 


9    ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

The  problem  of  estimating  the  log-survivor  function  for  a  fixed  time  t  of  a  first 
passage  time  for  a  semi-Markov  process  is  considered.  A  nonparametric  estimator  based  on 
an  exponential  approximation  to  the  survivor  function  is  introduced  and  its  asymptotic 
properties  are  studied.  Small  sample  behavior  of  this  estimator  is  studied  using  simulation 


20    DISTRIBUTION    AVAILABILITY  OF  ABSTRACT 
LXl  UNCLASSIFIED. UNLiM^ED       □  SAME  AS  RPT 


□  DTIC  USERS 


21     ABSTRACT  SECURITY   CLASSIFICATION 

UNCLASSIFIED 


?2a    NAME  OF  RESPONSIBLE  INDIVIDUAL 

P.    A.    Jacobs 


22b  TELEPHONE  (Include  Area  Code) 


(408)646-2258 


12<     OFFICE   SYMBOl 

55Jc 


)DFORM  1473.  84  mar 


S3  APR  ed  ton  may  be  used  until  exhausted 
a  i  ither  editions  are  obsolete 


SECURITY  CLASSIFICATION  OF   THIS   I 


ft  U  S    Govern. ni  nt   Prtntinq  Office     1986     606243 


18.      bootstrap;   jackknife 


1.   INTRODUCTION 

Finite  state  space  semi-Markov  processes  find  application  in  many  areas 
including  queueing,  reliability,  and  clinical  trials;  (cf  Cinlar  [1975],  Weiss 

and  Zelen  [1965],  Cox  [1985]).  Often  interest  in  the  application  centers  on 
the  distribution  of  a  first -passage  time  to  a  state  or  a  set  of  states 
representing  for  example  the  lifetime  of  a  system  or  the  end  of  a  busy  period 
of  a  server. 

In  this  paper  we  consider  the  problem  of  the  estimation  of  the  log 
survivor  function  of  a  first  passage  time  for  a  fixed  time  t.  We  suppose  that 
observational  data  are  known  about  the  semi  Markov  process  in  question. 
Although  the  approach  and  results  are  given  concretely  for  one  particular 
semi -Markov  process,  they  apply  more  widely. 

To  be  specific,  consider  N  individuals.  Let  X  (i)  be  the  state  of  ith 
individual  at  time  t.  We  will  assume  {X  (i);t>0}  are  independent  identically 
distributed  processes  having  the  same  probability  law  as  {X  ;t>0}.  The 
process  {X  ;t>0)  is  a  semi-Markov  process  with  three  states  {0,1,2}.  State  0 
is  absorbing  and  X  =1  .  The  sojourn  time  in  state  i  has  a  distribution 
function  F..  Upon  leaving  state  1,  the  process  transitions  to  state  0  with 
probability  9  and  to  state  2  with  probability  1-0.  From  state  2  the  process 
transitions  to  state  1  with  probability  1. 


Let 


D  =  inf  {t>0:   X   -  0} 


the  entrance  time  to  state  0 


The  problem  is  to  estimate  the  logarithm  of  P{D>t},  for  fixed  t,  from  data 
obtained  by  observing  the  N  individuals.  Three  possible  estimation  procedures 
will  be  considered.  The  three  estimators  use  different  amounts  of  information 
concerning  the  process.  One  procedure  uses  only  the  observed  entrance  time  to 
state  0  for  the  N  individuals  to  estimate  the  empirical  distribution  function 
of  D.  Another  procedure  makes  parametric  assumptions  concerning  the  sojourn 
time  distributions  F  and  uses  maximum  likelihood  to  estimate  the  probability. 
A  third  approach  uses  an  exponential  approximation  to  P{D>t}  and  empirical 
distributions  to  estimate  F  . 

In  Section  2,  the  estimation  procedure  based  on  an  exponential 
approximation  to  P{D>t}  is  described  and  asymptotic  results  are  obtained  for 
it. 

In  Section  3  the  other  two  estimation  procedures  are  described  and  the 
results  of  a  simulation  experiment  are  given.  The  experiment  was  conducted  to 
compare  the  performance  of  confidence  intervals  for  the  three  estimators  for 
small  and  moderate  numbers  of  individuals.  It  is  found  that  the  procedure 
based  on  the  maximum  likelihood  is  the  most  efficient  if  the  underlying  model 
is  correct.  However,  maximum  likelihood  confidence  interval  coverage  is 
sensitive  to  incorrect  model  assumptions.  The  nonparametric  confidence 
intervals  based  on  the  exponential  approximation  have  the  correct  coverage  as 
long  as  the  time  t  is  not  too  small.  The  conservative  binomial  confidence 
intervals  based  only  on  the  observed  times  of  absorption  are  the  least 
efficient  and  tend  to  overcover. 

2.   AN  ASYMPTOTIC  RENEWAL  ESTIMATE 

The  probability  P{D>t}  for  the  model  described  in  the  Introduction 


satisfies  the  equation 


P{D>t}  =  g(t)  +  (1-8)   /o  (Fj*F2)(ds)  P<D>t.-s}  (2.1) 


with 


g(t)  =  Fj(t)  +  (1-9)  /0  F^ds)  F2(t-s).  (2.2) 


This  is  a  renewal  equation  with  defective  inter-renewal  distribution 


L(t)  -  (l-e)(F1*F2)(t).  (2.3) 


Following  Feller  [1966]  page  326,  assume  there  exists  k    which  satisfies 
the  following  equation 


/°°  eKt   L(dt)  =  (1-9)  ♦  (k)  ♦  .(«)  =  1  (2.4) 

0  l  d 


where 


♦  .(C)  =  J°°  eCt  F  (dt).  (2.5) 

1       0 


If  (F  *F  )  is  not  arithmetic,  then  under  certain  integrabil i ty  conditions 

lim  eKt  P{D>t)  =  -  (2.6) 

t-^» 


where 


fJ   =   I°°     se"S    L(ds)  (2.7) 


and 


b  =  /°°  eKS  g(s)ds  =  %  *(k).  (2.8) 

0  K      * 


Following  an  approach  that  was  found  useful  in  another  context,  (Gaver 
and  Jacobs   ([1986]),   an  estimator  for  P{D>t}  can  be  obtained  by   using 

estimates  of  k  ,  b  and  ^    in  (2.6).   Specifically,  let  (S  (i),n=l M.}  be 

the  collection  of  sojourn  times  in  state  i  for  the  N  individuals  and  let  R  be 
the  number  of  transitions  from  state  1  to  state  2  for  all  the  N  individuals. 


N 
Put  8  =  —   and 

i  (O  -   1   z1  e?Sk(i).  (2.9) 

i   k=l 

i  =  1,2. 


An  estimator  of  k    is  the  solution,  k ,    to  the  equation 


(1-9)   ^(O^U)  =  1.  (2.10) 


Notice  that  the  left  hand  side  of  equation  (3.10)  is  nondecreasing  in  (  and  at 
C  =  0  equals  (1-9)  <  1.  Hence  (2.10)  has  a  unique  solution  which  can  be  found 
numerically . 


The  estimators  for  the  other  quantities  in  (2.6)  are 


.   ,      Ml     M2 
^  =  (1"e)M~M~   Z      Z      exP{«(sk(D+Sj(2))}  (SR(1)  +  sj(2))   (2.11) 
1  2   k=l    j-1 


b  =  t  ♦  (K).  (2.12) 

K 


The  asymptotic  exponential  estimator  for  P{D>t}  is 


PA(D>t}  =  t  e'Kt.  (2.13) 


2 . 1   Large  Sample  Properties  of  *  and  (n   P  (D>tl 


Let 


f(C.a)  =  1    (l-a)*1(e)*2(C)  =   1  ~(l-a)h(C) 


Expand  h(x)  in  Taylor  series  about  the  solution  k   of  (2.4).   Since  f(/c,0)  =  0, 


0  =  f{K,9)    =    1    (1-e)  [h(«)  +  (K-K)h'(K)    +  Ukx)2   h"(a«)]     (2.14) 


f(*,e)  +  (o-e)h(K)       (i-B)(K-K)[h'  (k)   +  i(x-/c)h"  (a*)] 


Thus 

f(«,e)  t  (e-8)  h(/<) 

K  K  =  (2.15) 

B 


where 


B  =  (1  e)[h'  («)  +  2(K  K)    h'  '  (**)] 


(2.16) 


and  a  is  a  constant  and  9    is  the  probability  of  a  transition  from  state  1  to 
state  0.   By  the  strong  law  of  large  numbers,  as  N~««>  f(*,0)-»O  and  9~9-*0. 
Further,  as  N-*»,  with  probability  one 


and 


h'(«)  -  (I-*)  J~  se  "S(F  *F  )(ds) 

0         l      6 


h(K)  -  ♦^k)  *2(k) 


(2.17) 


(2.18) 


Therefore,  «-«-»()  as  N-»«>,  with  probability  1  so  the  estimator  k    is  consistent. 

To  obtain  a  central  limit  theorem  for  «,  the  following  moment  condition 
is  needed.   Assume 


I      e    F  *F  (ds)  <  «. 

0      *  d 


(2.19) 


Since  f(*,9)    =  0 


y~F7  [  f(*,«)]  =  /T  [  f(«,e)   f(«,0)] 


N  (1-9)  [♦1(k)*2(k)   ♦1(«)*2(ic)] 


N  (1-0)  {[♦1(k)  -  ♦j(k)]  *2(k)  +  ♦1(k)[*2(#c)^2(k)]}.   (2.20) 


Thus,   (2.19)  and  the  central   limit  theorem  imply  that  /~N  [f(*,0)]   is 
asymptotically  normal  as  N-»°°  with  mean  0  and  variance 


K  T  K  T 

Varf  -  (1-9)    9{{l-t)*2(K)2   Var[e   !]  +  ^(k)2  Var[e   2]} 


(2.21) 


where  T.  is  a  random  variable  with  distribution  F  . 

i  l 

Further   the   central   limit   theorem   implies   that   as   N-»»,  J~H     [9-9]      is 

2      3 
asymptotically  normal  with  mean  0  and  variance  V.   =  9    (1-9)       +■  0  (1-0). 

b 

Hence,  as  N-*»,  Vn  (kk)  is  asymptotically  normal  with  mean  0  and  variance 


k(T,  *T, 


V«  =  (Vf  *  (V°  *2(/<))   Vb,[{1"9)  E[(T^T2'e    1   2  ]]        (2.22) 


A  similar  argument  shows  that  £ n  P  {D>t}  is  also  asymptotically  normally 
distributed  as  N-^>. 


3.   SIMULATION  RESULTS 

A  simulation  experiment  was  carried  out  to  investigate  small  and  moderate 
sample  size  behavior  of  confidence  interval  procedures  for  three  estimators 
of  (n    P(D>t}.    The  following  three  estimators  use  different  amounts  of 

information  about  the  semi-Markov  processes  {X  (i);  t>0},  i=l N,  for  the 

N  individuals 

3.1   The  Binomial  Estimator 

An  estimator  which  uses  only  the  observed  entrance  times  into  state  0  for 
the  N  individuals,  d  d    . .,d   is 


N 

^n  <[,>t>  "I  £   1<t.-)(V  <3 ■" 

1=1 


where 


1(t.«)(di) 


1    if  d.>t, 

1 


0   otherwise 


The  distribution  of   the  estimator  P  {D>t}   is   binomial   and  a  binomial 

D 

confidence  interval  [  L  ,  U  j  can  be  constructed  for  P{D>t}  cf.  Larson  [1982]  p 

B   B 

397.   The  binomial  confidence  interval  for  *n  P{D>t}  is  [( n   LD,fn  U  1  . 

D        D 


3.2  The  Asympiotic  Exponential  Estimator. 

The  asymptotic  exponential  estimator  for  P{D>t}  is  given  by  P  {D>t}  of 

(2.13).   The  estimator  for  the  en    P{D>t}  is  (r\    PA{D>t)  and  is  asymptotically 

A 

normal  as  N-»«>.  However,  since  the  simulation  experiment  is  for  small  to 
moderate  numbers  of  individuals,  two  nonparametric  methods  are  used  to  obtain 
confidence  intervals. 

The  jackknife  is  a  procedure  originally  introduced  by  Quenouille  [1956] 
for  bias  reduction,  and  adapted  by  Tukey  [1958]  to  obtain  approximate 
confidence  intervals.  Suppose  interest  is  on  a  parameter  p  that  is  estimated 
by  p  using  a  complex  calculation  from  data.  The  idea  is  that  of  assessing 
variability  by  recomputing  /3  after  removing  independent  subgroups  of  data  and 
then  using  the  recomputed  p  values  to  estimate  a  variance  which  is  in  turn 
applied  to  state  a  two-sided  confidence  interval  that  contains  the  true  p   with 


specified  confidence;  see  Efron  [1982 J  and  his  more  recent  work,  or  Mosteller 

and   Tukey   [1977]   for   more   details.     For   simulations   involving   N=20 

individuals,  the  jth  subgroup  consists  of  all  data  corresponding  to  the  jth 

individual.    For  simulations  involving  N=50  individuals  the  first  subgroup 

consists  of  all  data  corresponding  to  the  first  five  individuals  the  second 

subgroup  of  all  data  corresponding  to  the  second  five  individuals,  etc.   Some 

simulations  for  N=50  individuals  were  run  with  the  jth  subgroup  consisting  of 

the  data  for  the  jth  individual  alone;  according  to  some  work  these  should  be 

an  improvement  on  those  with  5  in  a  subgroup.    However,   the  resulting 

confidence  intervals  differed  little  from  those  obtained  by  leaving  out  5 

individuals  at.  a  time. 

The  bootstrap  is  an  alternate  method  that  may  be  used  to  assess  uncertainty 

of  estimators;  cf.  Efron  [1982]  and  his  more  recent  work.   In  the  simulations 

reported  here  a  bootstrap  replication  is  generated  as  follows.   Let  R  be  the 

n 

number  of  transitions  from  state  1  to  state  2  for  individual  n.   Let  {S  (i)} 

K 

be  the  collection  of  all  sojourn  times  in  state  i  for  all  individuals.   To 

generate  bootstrap  data  for  one  individual,  one  observation  is  drawn  at  random 

with  replacement  from  {R  };  call  the  observation  r;  r  +  1  observations  are  drawn 

n 

at  random  with  replacement  from  the  collection  of  state-1  sojourn  times 
{S.  (1)};  r  observations  are  drawn  at  random  with  replacement  from  the 
collection  of  state  2  sojourn  times  {S  (2)}.   Bootstrap  data  is  generated  for 

K 

N  individuals  and  the  estimate  tx\  P  {D>t}  is  computed  using  this  data.  This 
completes  one  bootstrap  replication.  B  =  100  bootstrap  replications  are  done. 
The  B  =  100  bootstrap  estimates  of  *n  P  {D>t}  are  ordered.  A  100  x  (lot)  % 
confidence  interval  is  constructed  using  the  a/2  and  l-a/2  quantiles  of  the 
bootstrap  estimates. 


3.3  The  Maximum  Likelihood  Estimator 

In  what  follows,  we  assume  the  sojourn  time  in  state  i  has  an  exponential 

distribution   with   mean   — .     The   log  likelihood   function   under   these 

i 


assumptions  is 


Z   [R  (n    (1-9)    +   en  9    «•  R  *  n  p   ♦  (1+R  )*n  p 

II  \  t  c.  II  J 

n=l 


(3.2) 


Vn(1)   Vn(2)] 


where  R   is  the  number  of  transitions  from  state  1  to  state  2  for  individual  n 
n 

and  T  (i)  is  the  total  time  spent  in  state  i  by  individual  n.   The  maximum 


likelihood  estimators  are 


N 


N  +  R 


(3.3) 


N+R 
*1  :  T(l) 


(3.4) 


R 


P2  :  TT2) 


(3.5) 


N 


N 


where  R  =  Z   R  ,  T(l)  =  Z  T  (1),  and  T(2)  =  Z  T  (2) 
n  n  n 


n=l 


n=l 


n=l 


Let  D  be  the  entrance  time  to  state  0  for  an  individual;  then 


A  +p  A  +p  9p 

S(t)  =  P{D>t}    {-J- 1  exp{A2t}   -L-l   exp{A1t}}  j—^ 


1 


1   2 


(3.6) 


10 


where  A,.A2  are  tne  roots  of  tne  equation 


9p1p2  f  y(Pj+P2)  +  y    =  0-  (3.7) 


The  maximum  likelihood  estimator  for  en  P{D>t}  can  be  obtained  from  (3.6)  and 
(3.7)  by  replacing  p  ,  p  and  9  by  their  maximum  likelihood  estimators  and 
taking  the  logarithm  of  the  resulting  estimator  of  F'{D>t}.  The  distribution 
of  the  resulting  estimator  is  asymptotically  normal  as  N-*».  Since  the  maximum 
likelihood  estimators  9  ,p  ,  and  p  are  uncorrelated ,  the  asymptotic  variance 
is 


2       ^       2       *  2 

Var[*n  P{D>t}]  -    *    [Var(3)  (||)   +  Var(pW^L)   *  Var(p  W^-)  ].  (3.8) 

ae         i  ap  2  ap 

sur  ^ 


Confidence  intervals  are  constructed  using  the  asymptotic  normality  of  the 
maximum  likelihood  estimator  and  the  asymptotic  variance  with  parameter  values 
being  replaced  by  their  corresponding  estimates.  Confidence  intervals  based 
on  the  maximum  likelihood  ratio  test  may  also  be  constructed  [cf.  Cox  and 
Hinkley  (1974)  pp.  343].  While  these  intervals  are  asymptotically  equivalent 
to  the  normal  confidence  intervals  that  are  computed,  they  may  be  better  for 
sample  sizes. 

3 .4  The  Simulation  Experiment 

All  simulations  are  carried  out  on  an  IBM  3033AP  computer  at  the  Naval 
Postgraduate  School  using  the  LLRANDOMII  random  number  generating  package;  [cf 


1  1 


Lewis  and  Uribe  (1981)  J.  Some  details  of  the  simulation  are  given  below.  A 
more  complete  account  can  be  found  in  Kim  [1987].  In  each  replication 
estimates  and  confidence  intervals  are  computed  for  the  logarithm  of  P { I)>  t } 
using  the  procedures  described  in  subsections  3.1-3.3.  Each  simulation  is 
replicated  300  times.  For  each  procedure,  the  number  of  confidence  intervals 
covering  the  true  value  of  Cn  P{D>t}  is  recorded.  Also  recorded  are  the 
number  of  intervals  that  are  too  low  (true  (n  P{D>t}  >  U,  the  upper  endpoint 
of  the  interval)  and  too  high  (true  (n  P(D>t}  <  L,  the  lower  endpoint  of  the 
interval).   The  average  length  of  the  confidence  intervals  is  also  computed. 

Tables  1-8  present  results  from  simulations  for  N=20  and  N=50 
individuals.  The  first  row  in  each  table  lists  the  confidence  interval 
procedures:  the  binomial  confidence  interval  for  logarithm  of  estimate  (3.1) 
(BIN);  the  maximum  likelihood  interval  (MLE);  the  jackknife  confidence 
interval  for  the  asymptotic  estimate,  (JK);  and  the  bootstrap  confidence 
interval  for  the  asymptotic  estimate,  (BT).  Nominal  80%  and  90%  confidence 
interval  for  (n  P{D>t}  are  computed  using  each  procedure  for  each  time 
t  =  0.5,  1.0,  1.5,  2,  3,  4.  Confidence  limits  that  are  greater  than  0  are  set 
equal  to  0. 

Tables  1-4  present  results  for  the  case  in  which  the  sojourn  times  in 
state  1  have  an  exponential  distribution  with  p  =1  and  the  sojourn  times  in 
state  2  have  an  exponential  distribution  with  mean  p  =0.1.  The  values  of 
P(D>t}  for  the  model  for  t  =  0.5,  1.0,  1.5,  2,  3,  4  are  given  in  Appendix  A. 

Tables  1-2  show  results  for  simulations  with  N=20  individuals. 
Table  1  shows  coverage  results.   Under  each  procedure  is  given  the  number  of 
intervals  that  cover  the  true  value  (C),  the  number  of  intervals  that  are  too 
low  (L)  and  the  number  of  intervals  that  are  too  high  (H).   In  parentheses 
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next  to  these  numbers  are  given  the  corresponding  fractions.  If  the  true 
confidence  level  of  a  confidence  interval  procedure  is  80%,  (respectively 
90%),  then  a  95%  confidence  interval  for  the  fraction  of  the  300  replications 


that  cover  the  true  value  is  0.80  ±  1.96     — —  (.2)  (.8)  =  [755,  .845], 

\j\J\J 


J — I 

(respectively  .9  ±  1.96   -^  (.1)(.9)  =  [.866,  .934]). 

Table  2  shows  the  average  length  of  confidence  intervals  for  £n  P{D>t) 
for  each  t  and  each  procedure.  In  parentheses  under  the  average  length  is  the 
estimated  standard  error  of  the  average  length. 

The  conservative  binomial  procedure  tends  to  overcover  and  has  the 
largest  average  confidence  interval  length.  The  maximum  likelihood  procedure 
has  the  correct  coverage  and  smallest  average  length.  The  confidence 
intervals  based  on  the  asymptotic  renewal  estimate  show  the  correct  coverage 
for  the  times  larger  than  t  =  0.5  and  1.0.  When  the  jackknife  and  bootstrap 
intervals  do  not  cover  they  tend  to  be  too  high.  The  average  length  of  the 
jackknife  and  bootstrap  intervals  is  between  those  of  the  binomial  and  maximum 
likelihood  intervals.  The  average  length  of  the  jackknife  intervals  is 
greater  than  those  for  the  bootstrap,  suggesting  that  the  latter  is  preferred; 
however  computational  cost  may  be  higher. 

Tables  3  and  4  present  coverage  results  and  average  confidence  interval 
lengths  for  a  simulation  experiment  with  N=50  individuals.  The  coverage 
results  are  similar  to  those  of  Table  1.  As  anticipated,  the  average 
confidence  interval  lengths  tend  to  be  smaller  for  N  =  50  individuals  than  for 
N=20  individuals.  The  average  bootstrap  interval  length  is  closer  to  that  of 
the  MLE  interval  for  times  larger  than  1  for  N=50  individuals  than  for  N^20 
individuals . 

Tables  5-8  report  results  of  a  simulation  experiment  in  which  the  sojourn 

1) 


times  in  states  1  and  2  have  the  following  distributions.  The  sojourn  time  in 
state  2  has  an  exponential  distribution  with  mean  —  =  0  1  as  before.   The 

sojourn  time  in  state  1  is  the  sum  of  two  independent  exponential  random 
variables  each  having  mean  -^.  Thus,  the  mean  sojourn  time  in  each  state  is 
the  same  as  for  the  model  used  in  Tables  14.  Values  of  P{D>1}  for  various 
values  of  t  for  this  model  are  given  in  Appendix  B.  Estimators  and  confidence 
intervals  for  fn  P{D>  }  are  computed  as  before.  In  particular-  the  maximum 
likelihood  estimators  assume  that  the  sojourn  times  in  state  1  have  an 
exponential  distribution  with  mean  1.  The  number  of  individuals  is  N=20  for 
the  simulation  in  Tables  5-6.  There  are  N=50  individuals  in  the  simulations 
of  Tables  7-8.  The  coverage  results  contained  in  Tables  5  and  7  show  that  the 
maximum  likelihood  confidence  intervals  are  sensitive  to  underlying  model 
assumptions.  Incorrect  model  assumptions  lead  to  both  overcoverage  and 
undercoverage .  The  confidence  intervals  based  on  the  exponential 
approximation  have  the  correct  coverage  for  the  times  larger  than  0.5.  The 
average  length  of  the  jackknife  and  bootstrap  intervals  is  smaller  than  those 
for-  the  binomial.  The  average  length  of  the  bootstrap  interval  is  smaller 
than  those  for  the  jackknife. 

4.   CONCLUSIONS 

The  simulation  results  indicate  that  the  maximum  likelihood  estimator  is 
most  efficient  when  the  correct  model  is  used.  As  has  been  noted  elsewhere, 
the  MLE  estimator  is  rather  sensitive  to  incorrect  model  assumptions.  The 
simulation  results  for  the  asymptotic  exponential  estimator  are  somewhat 
surprising  in  that  it  appears  that  the  time  t  does  not  have  to  be  very  large 
in  order  for  the  asymptotic  exponential  estimator  to  work  well.   Both  the 
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nonparametric  bootstrap  and  jackknife  confidence  intervals  for  the  asymptotic 
exponential  estimator  have  the  correct  coverage  for  moderately  large  t.  As 
implemented  the  bootstrap  intervals  tend  to  have  a  smaller  average  length  than 
the  jackknife.  However,  neither  a  complete  jackknife  nor  a  bootstrap  with 
more  than  100  replications  has  been  done  due  to  the  massive  amount  of 
computation  involved.  As  expected,  the  binomial  estimator  produces 
conservative  confidence  intervals  which  tend  to  overcover.  However,  the 
binomial  confidence  intervals  are  appropriate  for  all  times  and  all 
semi-Markov  models. 
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Appendix  A 


Values  of  P{D>t}  for  the  semi-Markov  Model 
with  exponential  sojourn  time  in  state  1 


P^l. 


P2=10,  6=0.5 


time  t 

P{D>t> 

0.5 

.7866 

1.0 

.6203 

1.5 

.4891 

2.0 

.3857 

3.0 

.2399 

4.0 

.1492 
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Appendix  B 


Values  of  P{D>t}  for  the  semi-Markov  Model 
with  nonexponential  sojourn  time  in  state  1 


Pj=2. 


P2=10. 


P3=2,  9=0.5 


time  t 

P(D>t} 

0.5 

.8652 

1.0 

.6743 

1.5 

.5157 

2.0 

.3930 

3.0 

.2278 

4.0 

.1320 
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Table  1 

Coverage  Results  for  Confidence  Intervals  for  8r\   P{D>t} 
Exponential  Sojourn  Times 
N=20.  6=O.E>,    p  =1.0.  p2=10.0 


Time      Confidence  Coverage 
Level 


BIN 


MLE 


JK 


BT 


0.5 

80% 

H 

14( 

.05] 

24  ( 

.08] 

106{ 

.35] 

140( 

•47) 

C 

272  ( 

.91] 

251( 

.84] 

174( 

.58] 

160( 

.53) 

I. 

14( 

.05] 

25( 

.08] 

20  ( 

.07] 

0( 

.0) 

90% 

H 

4( 

.01] 

14( 

.05] 

78  ( 

.26] 

99  ( 

.33) 

C 

282  ( 

.94] 

276  ( 

.92] 

211{ 

.70] 

201  ( 

.67) 

L 

14( 

.05] 

10{ 

.03] 

1    H< 

.04] 

0( 

.00) 

1.0 

80% 

H 

17( 

.061 

24  ( 

.081 

72  ( 

.241 

66  ( 

.22) 

C 

250  ( 

.83] 

1   250( 

.831 

210| 

'.701 

I   224  ( 

.75) 

1. 

33  ( 

.11] 

26  < 

;.09] 

)    18| 

[•06] 

10| 

[03) 

90% 

H 

7( 

.021 

1    14| 

.051 

1    46| 

'.15' 

►   33| 

[11) 

C 

275  ( 

.921 

|  277| 

.921 

1   245| 

r.821 

►   254| 

[•86) 

L 

18( 

.06; 

)    9| 

[.03] 

1     91 

[03; 

1    8| 

;.03 

1.5 

80% 

H 

15( 

.05' 

I   24< 

r.08' 

>    48| 

[.16' 

)    40 

[.13) 

C 

265  ( 

.88 

>  250) 

r.83 

)   2471 

r.82 

)   236 

[.79) 

L 

20  ( 

.07 

>   26| 

[.09 

)     5 

[.02 

)   24 

[.08) 

90% 

H 

15( 

.05 

>   14 

[.05 

)    33 

[.11 

)    19 

[.06) 

C 

278  ( 

.93 

)     277 

[.92 

)   266 

[.89 

)   266 

[.89) 

L 

7( 

.02 

)    9 

[.03 

)     1 

[.00 

)    15 

(.05) 

2.0 

80% 

H 

H( 

.04 

)   24 

[.08 

)    50 

[.17 

)   27 

(.09) 

C 

267  ( 

.89 

)  250 

[.83 

)   245 

[.83 

)   240 

(.80) 

L 

22  ( 

.07 

)   26 

[.09 

)     5 

[.02 

)   33 

(.11) 

90% 

H 

H( 

.04 

)   H 

[.05 

)    28 

[.09 

)   io 

(.03) 

C 

280  ( 

.93 

)  277 

[.92 

)   271 

[.90 

)   267 

(.89) 

L 

9( 

.03 

)    9 

[.03 

)     1 

[.00 

)   23 

(.08) 

3.0 

80% 

H 

23( 

.08 

)   23 

[.08 

)    40 

[.13 

)   19 

(.06) 

C 

271  ( 

.90 

)   251 

[.84 

)   252 

[.84 

)   238 

(.79) 

L 

6( 

.02 

)   26 

[.09 

)     8 

[.03 

)   43 

(.14) 

90% 

H 

8( 

.03 

)    14 

[.05 

)    22 

[.07 

)    9 

(.03) 

C 

286  ( 

.95 

)   277 

[.92 

)   276 

[.92 

)   260 

(.87) 

L, 

6( 

.02 

)    9 

[.03 

)     2 

[.01 

)   31 

(.10) 

4.0 

80% 

H 

19( 

.06 

)   23 

[.08 

)    36 

[.12 

)    15 

(.05) 

C 

271  ( 

.90 

)  251 

[.84 

)   253 

[.84 

)   231 

(.77) 

L 

10( 

.03 

>   26 

[.09 

)    11 

[.04 

)   54 

(.18) 

90% 

H 

4( 

.01 

)    14 

[.05 

)    21 

[.07 

)    7 

(.03) 

C 

289  ( 

.96 

)   277 

[.92 

)   277 

[.92 

)   261 

(.87) 

L 

7( 

.02 

)    9 

[.03 

)     2 

(.01 

)   32 

(.11) 
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Table  2 

Average  Length  of  Confidence  Intervals  for  €n   P{D>t} 
Exponential  Sojourn  Times 
N=20.  9=0.5.  p1=1.0,  p2=10.0 


Time      Confidence  BIN 
Level 


MLE 


JK 


BT 


0.5 

80% 

.370 

.144 

.208 

.185 

(.005) 

(.002) 

(.005) 

(.003) 

90% 

.460 

.186 

.271 

.245 

(.007) 

(.002) 

(.006) 

(.005) 

1.0 

80% 

.545 

.288 

.307 

.290 

(.008) 

(.004) 

(.006) 

(.005) 

90% 

.680 

.370 

.400 

.382 

(.009) 

(.005) 

(.008) 

(.006) 

1.5 

80% 

.716 

.431 

.467 

.436 

(.010) 

(.006) 

(.011) 

(.008) 

90% 

.893 

.553 

.608 

.579 

(.012) 

(.007) 

(.014) 

(.010) 

2.0 

80% 

.915 

.574 

.645 

.597 

(.016) 

(.008) 

(.015) 

(.011) 

90% 

1.14 

.737 

.840 

.793 

(.020) 

(.010) 

(.019) 

(.015) 

3.0 

80% 

1.38 

.860 

1.02 

.932 

(.033) 

(.011) 

(.024) 

(.018) 

90% 

1.71 

1.10 

1.33 

1.24 

(.038) 

(0.15) 

(.032) 

(.025) 

4.0 

80% 

2.11 

1.15 

1.40 

1.27 

(.069) 

(.015) 

(.033) 

(.026) 

90% 

2.58 

1.47 

1.82 

1.70 

(.076) 

(.020) 

(.044) 

(.034) 
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Table  3 


Coverage  Results  for  Confidence  Intervals  for  €n   P{D>t} 
Exponential  Sojourn  Times 
N=50.  6=0.5,  p1=1.0.  p2=10.0 


Time      Confidence  Coverage 
Level 


BIN 


MLE 


JK 


BT 


0.5 

80% 

H 

17(.06) 

27  ( 

.09 

)   H9( 

.40) 

1511 

[•50) 

C 

260(.87) 

245  ( 

.83 

)   167( 

.56) 

149| 

[•50) 

L 

23(.08) 

28( 

.09 

)    14( 

.05) 

0| 

[00) 

90% 

H 

5(.02) 

14( 

.05 

)   83( 

.28) 

110| 

[•37) 

C 

288(.96) 

274  ( 

.91 

)   209  ( 

.70) 

190| 

[•63) 

L 

7(.02) 

12( 

.04 

I    8( 

.03) 

o< 

[00) 

1.0 

80% 

H 

28(.09) 

27  ( 

.09 

)   70( 

.23) 

80| 

[•27) 

C 

245(.82) 

245  ( 

.82 

}  216( 

.72) 

213| 

[•71) 

L 

27(.09) 

28( 

.09; 

)    14( 

.05) 

7| 

[02) 

90% 

H 

8(.03) 

14( 

.05' 

\        45( 

.15) 

45| 

[15) 

C 

282(.94) 

274  ( 

.91' 

>   274  ( 

.82) 

250< 

[•83) 

L 

10(.03) 

12( 

.04; 

I    8( 

.03) 

5| 

[02) 

1.5 

80% 

H 

19(.06) 

27  ( 

.09' 

1   53( 

.17) 

46| 

[15) 

C 

262(.87) 

245  ( 

.821 

►   235( 

.78) 

233  < 

[•78) 

L 

19(.06) 

28( 

.09; 

1    12( 

.04) 

21| 

[07) 

90% 

H 

11(.04) 

14( 

.051 

)   32( 

.11) 

23| 

[08) 

C 

278(.93) 

274  ( 

.9r 

)     262  ( 

.87) 

266| 

[•89) 

L 

11(.04) 

12( 

.04; 

!    6( 

.02) 

HI 

[04) 

2.0 

80% 

H 

18(.06) 

27  ( 

.09' 

>   47  ( 

.16) 

31| 

[•10) 

C 

260(.87) 

245  ( 

.82' 

)   240  ( 

.80) 

243| 

[•81) 

L 

22(.07) 

28( 

.09; 

>   13( 

.04) 

26| 

[09) 

90% 

H 

9(.03) 

14( 

.05" 

!   22( 

.07) 

15| 

[05) 

C 

277(.92) 

274  ( 

.91' 

}  270  ( 

.90) 

2681 

[.89) 

L 

14(.05) 

12( 

.04; 

)    8( 

.03) 

17| 

[•06) 

3.0 

80% 

H 

20(.07) 

27  ( 

.09" 

)   36( 

.12) 

25| 

[08) 

C 

262(.87) 

245  ( 

.82" 

)  250( 

.83) 

2401 

[•80) 

L 

18(.06) 

28  ( 

.09; 

)   H( 

.05) 

35 

[12) 

90% 

H 

13(.04) 

14( 

.05' 

)   20( 

.07) 

HI 

[04) 

C 

283(.94) 

275  ( 

.92' 

1  274  ( 

.91) 

267| 

[•89) 

L 

4(.01) 

H( 

.04; 

I    6( 

.02) 

221 

[07) 

4.0 

80% 

H 

14(.05) 

27  ( 

.09' 

)   35( 

.12) 

20| 

[07) 

C 

278(.93) 

245  ( 

.82' 

|  248  ( 

.83) 

239| 

[•80) 

L 

8(.03) 

28( 

.09; 

1   17( 

.06) 

41| 

[•H) 

90% 

H 

6(.02) 

14( 

.05] 

19( 

.06) 

10| 

[03) 

C 

286(.95) 

275  ( 

.92] 

274  ( 

.91) 

262| 

[•87) 

L 

8(.03) 

H( 

.04] 

7( 

.02) 

28 1 

.09) 

2! 


Table  4 

Average  Length  of  Confidence  Intervals  for  £n   P{D>t} 
Exponential  Sojourn  Times 
N=50.  6=0.5.  p1=1.0.  p2=10.0 


Time      Confidence  BIN 
Level 


MLE 


JK 


BT 


0.5 

80% 

.216 

.089 

.155 

.127 

(.002) 

[.001)    | 

[.004) 

[.002) 

90% 

.271 

.114 

.205 

.162 

(.002)     I 

[.001)    | 

[.006) 

[.003) 

1.0 

80% 

.319 

.177 

.202 

.181 

(.003)     | 

[.001)    | 

;.004)    | 

[.002) 

90% 

.401 

.227 

.268 

.234 

(.003)     | 

[.002)    | 

[.005)    | 

[.002) 

1.5 

80% 

.416 

.245 

.282 

.262 

(.003)     < 

[.002)    | 

[.005)    I 

[.003) 

90% 

.524 

.340 

.374 

.340 

(.004)     { 

[.003)    | 

[.007)    | 

[004) 

2.0 

80% 

.524 

.353 

.377 

.352 

(.005)     | 

[.003)    ( 

[.007)    | 

[004) 

90% 

.628 

.453 

.500 

.458 

(.006)     | 

;.004)    ( 

[.009)    | 

[.005) 

3.0 

80% 

.750 

.529 

.585 

.545 

(.008)     ( 

[004)    ( 

;.on)       | 

[.006) 

90% 

.940 

.579 

.775 

.706 

(.009)     ( 

[.005)    | 

[•014)    1 

[.008) 

4.0 

80% 

1.05 

.705 

.799 

.738 

(.015)     ( 

.006)    ( 

[.015)    I 

[.009) 

90% 

1.31 

.905     1 

1.06 

.959 

(.019)     ( 

.007)    | 

[.020)    | 

[011) 

22 


Table  5 

Coverage  Results  for  'Confidence  Intervals  for  in   P{D>t} 

Nonexponential  Sojourn  Times 

N=20.  9=0.5.  p1=2.0.  p  =10.0.  p3=2.0 


Time      Confidence  Coverage 
Level 


BIN 


MLE 
(Exp) 


JK 


BT 


0.5 

80% 

H 

12( 

.04] 

0( 

.0) 

105( 

.35] 

135( 

.45) 

C 

278  ( 

.93] 

53  ( 

.18] 

183( 

.61] 

162( 

.54) 

L 

10( 

.03] 

1      247  ( 

.82] 

12( 

.04] 

3( 

.00) 

90% 

H 

0( 

.OO1 

0( 

'.001 

73  ( 

.24] 

82  ( 

.37) 

C 

290  ( 

.97] 

1      109( 

.36] 

220  ( 

.73] 

217( 

.72) 

L 

10( 

;.03] 

1      191{ 

[•64] 

7( 

.02] 

1( 

.0) 

1.0% 

80% 

H 

16| 

'.051 

1          K 

•0) 

45  ( 

.15' 

32  { 

•11) 

C 

254  j 

'.851 

)     221  ( 

.741 

1        240  ( 

.801 

1        241  ( 

.80) 

L 

30| 

[•10] 

1        78( 

[•26] 

1          15( 

.05] 

27 1 

[09) 

90% 

H 

5| 

[•02] 

>          0| 

[0) 

29 1 

.10] 

1          111 

[04) 

C 

285| 

'  .951 

|      266| 

'.891 

I        263  ( 

'.881 

1       271( 

[•90) 

L 

10| 

;.03] 

1        34  j 

[•11] 

1            8< 

[-02] 

1          18! 

[06) 

1.5 

80% 

H 

21| 

r.07' 

)          8! 

'.OS1 

I          40| 

'.131 

1          23) 

[08) 

C 

260| 

r.87' 

>     257< 

r.86' 

I        244( 

.811 

I        232( 

[.77) 

L 

19! 

[•06] 

1       35) 

[12] 

1          16| 

[05] 

I          45( 

[•15) 

90% 

H 

5| 

.021 

1         2| 

[.01] 

1          24( 

'.08' 

1            7| 

[02) 

C 

276) 

.921 

>     283< 

[.941 

I        273( 

'.911 

1        264| 

[•88) 

L 

19| 

[•06] 

1        15( 

[05] 

1            3| 

[01] 

I          29| 

[.10) 

2.0 

80% 

H 

111 

r.04" 

1        17< 

r.06' 

1          39! 

.13 

)           18| 

[06) 

C 

270| 

r.90' 

>     260| 

r.871 

>        246( 

.82' 

)        232| 

[.77) 

L 

19< 

[-06] 

)        23! 

[•08] 

1          15( 

[-05] 

)          50| 

[.17) 

90% 

H 

111 

'.04' 

1          7< 

r.02' 

1          19( 

r.06" 

>            7! 

[02) 

C 

284| 

r.95' 

|      284| 

r.95' 

)        277( 

r.92' 

}        263 

[.88) 

L 

5| 

[02] 

1          9| 

[-03] 

)            4! 

[01] 

)          30 

[10) 

3.0 

80% 

H 

13| 

\04' 

}        31| 

MO' 

)          37! 

.12 

)           17 

[06) 

C 

273! 

r.91" 

)     254 

r.85 

}        2511 

r.84 

)        230 

(.77) 

L 

14! 

[-05] 

I        15 

[.05 

}          121 

[.04 

)          53 

(18) 

90% 

H 

51 

r.02 

)        17 

[.06 

)          19 

[.06 

)            6 

(.02) 

C 

281 

r.94 

)     279 

[.93 

)        277 

[.92 

)        261 

(.87) 

L 

14 

[.05 

)           4 

[.01 

>            4 

[.01 

)          33 

(.11) 

4.0 

80% 

H 

4 

r.01 

)        44 

(.15 

)          38 

(.13 

)          1^ 

(.06) 

C 

279 

r.93 

)      245 

[.82 

)        251 

(.84 

)        229 

(.76) 

L 

17 

[.06 

)        11 

(.03 

)          11 

(.04 

)          54 

(.18) 

90% 

H 

4 

r.01 

)        23 

(.08 

)          19 

(.06 

)            6 

(.02) 

C 

296 

r.98 

)     273 

(.91 

)        277 

(.92 

)        260 

(.87) 

L 

0 

[.00 

)           4 

(.01 

)             4 

(.01 

)          34 

(.11) 
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Table  6 

Average  Length  of  Confidence  Intervals  for  in   P{D>t} 
Nonexponential  Sojourn  Times 
N=20.  0=0.5.  p1=2.0.  p2=10.0  p3=2.0 


Time      Confidence  BIN        MLE       JK         BT 
Level  (Exp) 


0.5         80%     .289       .144  .141  .130 

(.004)  (.002)  (.003)  (.003) 

90%     .359       .185  .183  .175 

(.005)  (.002)  (.004)  (.004) 

1.0         80%     .485       .286  .269  .253 

(.007)  (0Q3)  (.006)  (.004) 

90%     .604       .367  .351  .331 

(.008)  (.004)  (.008)  (.006) 

1.5         80%     .681       .428  .451  .419 

(.010)  (.005)  (.010)  (.007) 

90%     .850       .550  .588  .552 

(.013)  (.007)  (.014)  (.010) 

2.0         80%     .884       .571  .643  .594 

(.013)  (.007)  (.015)  (.017) 

90%    1.10        .733  .838  .786 

(.017)  (.009)  (.019)  (.015) 

3.0         80%    1.46        .855  1.03  .953 

(.034)  (.010)  (.024)  (.018) 

90%    1.81  1.10  1.35  1.26 

(.043)  (.013)  (.032)  (.025) 

4.0         80%    2.34  1.14  1.43  1.31 

(.079)  (.013)  (.034)  (.025) 

90%    2.84  1.46  1.86  1.74 

(.084)  (-017)  (.044)  (.035) 
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Table  7 

Coverage  Results  for  Confidence  Intervals  for  in   P{D>t} 

Nonexponential  Sojourn  Times 

N=50.  6=0.5.  p1=2.0.  p2=10.0,  p3=2.0 


Time      Confidence  Coverage    BIN 
Level 


MLE 
(Exp) 


JK 


BT 


0.5 

80% 

H 

18( 

.06] 

0(.0) 

130( 

.43] 

149( 

.50) 

C 

261  ( 

.87] 

3(.01] 

165( 

.55] 

150( 

.50) 

L 

21( 

.07] 

297(.99] 

5( 

.02] 

1( 

•0) 

90% 

H 

10( 

.03] 

0(.00] 

99( 

.33] 

106( 

.35) 

C 

283  ( 

.94] 

6(.02] 

200  ( 

.67] 

194( 

.65) 

L 

?( 

.02] 

294(.98] 

1( 

.00] 

0( 

■0) 

1.0 

80% 

H 

22  ( 

.071 

2(01] 

\          42( 

.141 

39  ( 

.13) 

C 

260  ( 

.871 

I   145(.48] 

|    242( 

'  .8V 

240  ( 

.80) 

L 

18( 

.06] 

I   153(.51] 

1    16| 

;.05] 

21( 

.07) 

90% 

H 

H( 

.041 

>   o(.oo; 

1    28| 

'  .091 

1    21( 

.07) 

C 

280  ( 

.93' 

|  213(.7T 

I   263| 

r.88' 

\       265  ( 

.88) 

L 

9( 

.03; 

)   87(.29; 

)     9| 

[03; 

1    14( 

.05) 

1.5 

80% 

H 

27  ( 

.09' 

)    6(.02 

)    311 

M0' 

)    25( 

.08) 

C 

254  ( 

.85' 

)  242(.81 

)   247| 

[.82 

>   243  ( 

.81) 

L 

19( 

.06 

)   52(.17 

}    22 

[.07 

>    32( 

.11) 

90% 

H 

H( 

.04 

)    4(.01 

)    21 

[.07 

>    19( 

.06) 

C 

276  ( 

.92 

)  276(.93 

)   271 

[.90 

)   262  ( 

.87) 

L 

13( 

.04 

)   20(.07 

)     8 

[.02 

)    19( 

.06) 

2.0 

80% 

H 

25( 

.08 

)   13(.04 

)    29 

[.10 

)    24( 

.08) 

C 

253( 

.84 

)  264( . 88 

)   250 

[.83 

)   239( 

.80) 

L 

22  ( 

.07 

)   23(.08 

)    21 

[.07 

)    37  ( 

.12) 

90% 

H 

16( 

.05 

)    7(.02 

)    20 

[.07 

)    12( 

.04) 

C 

274  ( 

.91 

)  282(.94 

)   269 

[.90 

)   264( 

.88) 

L 

10( 

.03 

)   11(.04 

)    11 

[.04 

)    24( 

.08) 

3.0 

80% 

H 

25( 

.08 

)   34(.ll 

)    31 

[.10 

)    21( 

.07) 

C 

247  ( 

.82 

)  253(.84 

)   248 

[.83 

)   241  ( 

.80) 

L 

28( 

.09 

)   13(.04 

)    21 

(.07 

)    38( 

.13) 

90% 

H 

17( 

.06 

}   20(.07 

)    18 

[.06 

)    H( 

.04) 

C 

273  ( 

.91 

)  275(.92 

)   271 

[.90 

)   266( 

.89) 

L 

10( 

.03 

)    5(.02 

)    11 

(.04 

)    23( 

.08) 

4.0 

80% 

H 

H( 

.04 

)   47( . 16 

)    31 

(.10 

)    20( 

.07) 

C 

266  ( 

.89 

)   246(.82 

)   244 

(.81 

)   238( 

.79) 

L 

23( 

.08 

)    7(.02 

)    25 

(.08 

)    42( 

.14) 

90% 

H 

6( 

.02 

)   29( . 10 

)    19 

(.06 

)    10( 

.03) 

C 

285( 

.95 

)  267(.89 

)   270 

(.90 

)   264( 

.88) 

L 

9( 

.03 

)    4(.01 

)    11 

(.04 

)    26( 

.09) 
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Table  8 


Average  Length  of  Confidence  Intervals  for  £n   P{D>t} 
Nonexponential  Sojourn  Times 


N=50.  6=0.5, 


P1=2.0. 


P2=10.0. 


p3=2.0 


Time      Confidence  BIN 
Level 


MLE 


JK 


BT 


0.5 

80% 

.169 

.089 

.091 

.081 

(.002)            | 

;.ooi)      \ 

\002)            ( 

.001) 

90% 

.211 

.114 

.121 

.104 

(.002)             | 

[001)          | 

;.003)            ( 

.002) 

1.0 

80% 

.284 

.176 

.167 

.153 

(.002)             | 

;.ooi)      | 

;.003)            ( 

.002) 

90% 

.357 

.227 

.221 

.199 

(.003)             | 

[.002)          | 

[004)            ( 

;.002) 

1.5 

80% 

.356 

.264 

.269 

.250 

(.003)            | 

;.002)          I 

[.005)            ( 

;.oo3) 

90% 

.498 

.339 

.357 

.324 

(.004)            | 

;.002)          I 

[.006)            I 

[004) 

2.0 

80% 

.514 

.352 

.378 

.351 

(.005)            | 

[.003) 

[.007)            | 

[004) 

90% 

.646 

.452 

.501 

.456 

(.006)            | 

[.003) 

[.010)            | 

[.005) 

3.0 

80% 

.784 

.527 

.600 

.558 

(.009)            | 

[•004) 

[.011)            | 

:.006) 

90% 

.983 

.677 

.795 

.726 

(.011)            | 

;.005) 

[.014)            | 

:.009) 

4.0 

80% 

1.16 

.703 

.824 

.765 

(.019)            | 

[.063)          | 

[.015)            | 

:.009) 

90% 

1.45 

.903 

1.09 

.998 

(.024)            ( 

[.064)          | 

[.020)            | 

[012) 
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